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ABSTRACT 



Aims. We have calculated variability Doppler boosting factors, Lorentz factors, and viewing angles for a large sample of sources by 
using total flux density observations at 22 and 37 GHz and VLBI data. 

Methods. We decomposed the flux curves into exponential flares and determined the variability brightness temperatures of the fastest 
flares. By assuming the same intrinsic brightness temperature for each source, we calculated the Doppler boosting factors for 87 
sources. In addition we used new apparent jet speed data to calculate the Lorentz factors and viewing angles for 67 sources. 
Results. We find that all quasars in our sample are Doppler-boosted and that the Doppler boosting factors of BL Lacertae objects are 
lower than of quasars. The new Lorentz factors are about twice as high as in earlier studies, which is mainly due to higher apparent 
speeds in our analyses. The jets of BL Lacertae objects are slower than of quasars. There are some extreme sources with very high 
derived Lorentz factors of the order of a hundred. These high Lorentz factors could be real. It is also possible that the sources exhibit 
such rapid flares that the fast variations have remained undetected in monitoring programmes, or else the sources have a complicated 
jet structure that is not amenable to our simple analysis. Almost all the sources are seen in a small viewing angle of less than 20 
degrees. Our results follow the predictions of basic unification schemes for AGN. 
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1. Introduction 

I All radio-bright active galactic nuclei (AGN) have relativistic 
■ jets emitting synchrotron radiation. The jets can at the sim- 
plest level be modelled by using two intrinsic parameters, the 
Lorentz factor (F), which describes the speed of the jet flow, 
' and the viewing angle (6), which is the angle between the jet 
I axis and the line of sight to the observer. These parameters can 
be calculated if the Doppler boosting factor (D) and the appar- 
ent speed ySapp - v/c are known. We can find out the /S^pp from 
Very Long Baseline Interferometry (VLBI) observations, and in 
the p ast few years major progr ess has been made in t his area 
(e.g. jjorstad et al . 2001; Hom an et all 120011; Ik ellerma iin et al.1 
l2004t iJorstad et al.ll2005i: IPiner et al.l l2007'; 'Bri tzen et al.ll2008l) . 
The Doppler boosting factors can be ca lculated in various ways, 
and d ifferent methods are compared in iLahteenmaki & Valtaoial 
([1999) (hereafter LV99). 

A common way to calculate the Doppler boosting fac- 
tors is to combine X-ray observations with VLBI compo- 



nent fluxes (e. g. iGhiselli niet al.1 [19931 iGuiiosa&Dalvl [19961 
iGuerra & Ddvlll99 7; Britze n et al.ll2007h . This method assumes 
inverse Compton (IC) origin of the X-ray emission, and that the 
same synchrotron photons forming the lower frequency radia- 
tion are responsible also for the IC emission. Assuming that the 
VLBI observations are done at the spectral turnover frequency, a 
predicted X-ray flux can be calculated. By comparing this to the 
observed X-ray flux, and by interpreting the excess flux as due 
to Doppler boosting, the Doppler boosting factors can be calcu- 
lated. If the VLBI frequency is not at the turnover, large errors 
are induced in the Doppler boosting factors. This method also 



suffers greatly from non-simultaneous X-ray and VLBI data, and 
as was argued in LV99, gives much less accurate estimates for 
the Doppler boosting factors. 

Using VLBI it is possible to directly observe the bright- 
ness temperature of the source (rb,obs)- This can be compared 
to the intrinsic brightness temperature of the source (Tb^int), 
which is often assumed to be the equipartiti on temperature 
(req) jReadheadl 119941: iLiihteenmaki et al.l [1999) . The excess 
of Throbs is interpreted as caused by Doppler boosting. This 
method also requires the values to be obtained at the turnover 
frequency, which enhances the errors in the Doppler boosting 
factors (LV99). 

Another way t o use VLBI observations is shown in 
iJorstad et al.| (l2005b who estimated the variability Doppler 
boosting factors of 15 AGN using Very Long Baseline Array 
(VLB A) data at 43 GHz. They calculated the flux decline time 
(Tots <^ Tj„,D) of a component in the jet and compared it to 
the measured size of the VLBI component (which does not de- 
pend on D). Assuming the intrinsic variability timescale corre- 
sponds to the light-travel time across the knot, they estimated the 
Doppler boosting factors. They also estimated the Lorentz fac- 
tors and the viewing angles for these sources by using apparent 
speed data. 

Variability timescales can also be obtained from total flux 
density (TFD) observations. This is the method used in LV99 
and we use the same method in our analyses. We decompose 
each flux curve into exponential flares and calculate the vari- 
ability timescale of each flare. From this we gain the observed 
brightness temperature, which is boosted by in comparison 
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with rb,int- This makes possible the extracti on of the variabil- 
ity D oppler factor Dyar if T\, int is known. In iLahteenmaki et all 
(11999 ) it was argued, based on observations, that in every large 
flare T\, m reaches the equipartition temperature Teq = 5 x lO'^/T. 

In LV99 a sample of 81 sources was studied at 22 and 
37 GHz frequencies. They calculated the Doppler boosting fac- 
tors based on observations from a period of over 15 years. For 
each source they combined the results of the two frequency 
bands and chose the fastest flare for the analysis. In addition they 
determined the Lorentz factors and the viewing angles for 45 
sources. We have done similar calculations for a larger sample 
of sources and using data from almost 30 years of monitoring. 
We will also compare our results with the results of LV99 to see 
if the values have changed during the past 10 years. 

The paper is organised as follows: in Sect. |2] we describe 
the source sample and the method used. In Sect. |3] we calculate 
the variability Doppler boosting factors and compare our results 
with LV99 and other related studies. The Lorentz factors and 
the viewing angles are presented in Sect. |4] and the discussion 
follows in Sect. |5] Finally conclusions are drawn in Sect. |6] 
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Fig. 1. Flux curve of the HPQ source 1 156-1-295 (points) decom- 
posed into exponential flares (solid line) at 22 GHz 



2. Data and the method 

Our sample consists of 87 bright, well-monitored AGN from 
the Metsahovi Radio Observatory monitoring list. This is a 
rather good approximation of a complete flux limited sample of 
the brightest compact northern sources. (Missing are a handful 
of sources for which we had insufficient data to calculate the 
Doppler factors.) The sources have been observed regularly for 
almost 30 years at 22 and 37 GHz with the Metsahovi 14m tele- 
scope (ISalonen et all 1 19871; iTerasranta et al.llT992l 1 19981 l2004l 
l2005h . Details of the observation method and data reduction pro- 
cess are described in ITerasranta et al.l (Il998h . Our study also in- 
cludes unpublished data at 37 GHz from December 2001 until 
the end of 2006. Data of BL Lacertae objects at 37 GHz from 
December 2001 until April 2005 are published in Nieppola et al.' 
12007). In our sample we have 30 high polarisation quasars 
(HPQs), which have optical polarisation exceeding 3 percent at 
some point in the past. In addition we have 22 low polarisa- 
tion quasars (LPQs), 8 quasars (QSOs) for which no polarisation 
data were available, 22 BL Lacertae objects (BLOs), and 5 ra- 
dio galaxies (GALs). We have used the classification used in the 
MOJAV^lpro gramme whenever possible. All the quasars in our 
sample can be classified as flat spectrum radio quasars (FSRQs), 
which show blazar-like properties. Our sample does not contain 
ordinary quasars that have larger viewing angles and which do 
not show the rapid variability common for FSRQs. In our anal- 
yses we have usually combined QSOs with LPQs, following the 
choice ofLV99. 

We decomposed the flux curves into exponential flares of the 
form 



AS(t) = 



ASr 
ASr 



,(f-W)/T^ 
C-'m«)/l'3T 



f < t„ 
t > t„ 



(1) 



where A5max is the maximum amplitude of the flare in janskys, 
/max is the epoch of the flare maximum and t is the rise time 
of the flare. The method is described in detail in IValtaoia et alj 
(I1999I) . An example of a fit is presented in Fig. [1] where the 
flux curve of a HPQ source 1156-1-295 is decomposed into ex- 
ponential flares. The solid line represents the fitted sum curve 
of the individual components and the points are actual observed 



data. We can see that the overal l correspondenc e betw een the 
two is very good. In LV99 and Savolainen et al. I (I2002h it was 
also shown that the individual exponential flare components cor- 
respond very well to the emergence of new VLBI components, 
indicating that these flares obtained from the fits are indeed re- 
lated to the actual jet physics. 

From the fits we obtain the necessary parameters to calcu- 
late the observed variability brightness temperature of the source 
T^b.var (in the source proper frame) 



rb,var= L548X 10" 



AS max"! 

v2t2(i +zy 



(2) 



where v is the observed frequency in GHz, z is the redshift, d]^ is 
the luminosity distance in metres, and AS' max and t are defined in 
Eq.[T] The numerical factor in Eq.|2]coiTesponds to using Hq - 
72kms"' Mpc"\ 0,„ = 0.27 and Qa - 0.73, and to assuming 
that the source is a homogeneous sphere. In the calculation of the 
luminosity distances we have made use of the python version of 
the co smology calculator created by Edward L. WrighQ ( Wrigh^ 
I2006I) . 

The variability Doppler factor can then be calculated as 



1/3 



(3) 



We used the value Tbant - 5 x 10 "^K. This value was 
suggested by Readhead ( 199 4j) an d the use of it was jus- 
tified in Lahteenmaki et all (Il999l) . Based on simulations, 
iKellermann et al. (2004F also find the Tb int to be of the order 
oflO"K. 

By combining the Doppler boosting factors with apparent su- 
perluminal component velocities ySapp, obtained using VLBI, we 
can calculate the variability Lorentz factors Fvai and the viewing 
angles ^vai by using Eqs.|4]and|5] 



•2 H- +\ 

app var 



2Dv 



arctan ■ 



2/3a 



app 



app 



(4) 
(5) 



' http://www.physics.purdue.edu/MOJAVE/ 



" http://www.astro.ucla.edu/ wright/CosmoCalc.html 
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We obtained 67 ySapp values from the MOJAVE sample, ob- 
served with the VLBA at 15 GHz. The values are taken from 
the website on September 9, 2008 and some of them may still 
be preliminary and change slightly in the final results (Lister et 
al. in preparation). This is the most homogeneous and largest 
sample of jSapp available at higher radio frequencies. These val- 
ues represent the fastest reliable speed in each jet as measured by 
the MOJAVE programme at 15 GHz, using VLBA data spanning 
between 5 to 13 years, depending on the individual source. 

3. The variability Doppler boosting factors 

3.1. Estimation of fhe Doppler boosting factors 

We were able to calculate the Doppler boosting factor (Dyar) for 
86 sources at 22 GHz, and for 72 sources at 37 GHz. For many 
sources we were able to determine the Dvai for more than one 
flare. The medians of Dvar are slightly larger at 22 GHz than at 
37 GHz, which could be due to different intrinsic brightness tem- 
peratures at the two frequency bands. In our analysis we chose 
the fastest flare of each source (either at 22 or 37 GHz, whichever 
had the fastest flare) to calculate the Dvai, as was done in LV99. 
The argument for using the fastest flare to determine Dvar is that 
they are most likely to reach the limiting brightness temperature 
and least likely to suffer from blending of flares which tends to 
increase the fitted timescale. This way we were able to calculate 
the Dysi for all the 87 sources in our sample. By visual examina- 
tion, we divided the fits into three categories based on the good- 
ness of the fit. No single numerical value, such as the test, is 
alone suitable for describing the goodness, because these usually 
characterise the entire flux curve, while we have only used one, 
fastest flare, to determine the ZJvai - In addition, relatively large 
error bars in some fainter sources cause the;^'^ value to be small, 
while we consider a fit to be better when the scatter among the 
datapoints is small. 

We classified the Dvai as excellent (E, 21 sources), good (G, 
24 sources) or acceptable (A, 42 sources). In the fits classified as 
excellent, the exponential decomposition follows the datapoints 
quite precisely, as in the flares after the year 1997 in Fig. [1] The 
fit is also unambiguous and other functions do not describe the 
behaviour as well. Larger flares in Fig. [T] before 1995 would 
mainly be defined as good, because in these flares the fit follows 
the flux curve well, but there is also some scatter, and in some 
cases there are not as many datapoints to define the fit as in the 
ones defined as excellent. In the fits classified as acceptable, the 
scatter around the fit is still larger. This is often the case when the 
flux level of the source is modest and the errorbars consequently 
large. However, none of our results change significantly if we 
exclude the acceptable sources. All Doppler boosting factors are 
shown in Table[T]where the B1950-name, other commonly used 
name, type of the object, frequency of the Dvar determination, 
redshift, quality of the Dvar, log Tb.var, £>var, jSapp, Fvar, Avar, core 
dominance parameter R, and maximum optical polarisation Pmax 
and its reference are listed. 

It is difficult to determine exact error estimates for the Dvar 
of each source. Some indications can be obtained from the stan- 
dard deviation of Dvar calculated from the various flares in one 
source. We calculated the deviations for all the sources classi- 
fied as excellent, which had more than one flare to determine 
the ZJvar, including also flares determined as good. There were 
45 such cases (including all the fits at both 22 and 37 GHz) and, 
on average, each source had 5.7 well-defined flares. The median 
standard deviation for these is ~ 27%, which can be thought of 
as an indication of the upper limit for the error estimate since in 



Table 2. Median values of log(rb,var) and D, 



Type 


N 


l0g(rb,va,)[/f] 




HPQ 


30 


14.31 


15.98 


LPQ 


30 


13.92 


11.90 


FSRQ 


60 


14.19 


14.61 


BLO 


22 


13.09 


6.25 


GAL 


5 


11.65 


2.07 


ALL 


87 


13.94 


12.02 



many cases the change in the Dyai of individual flares can also be 
due to differences in the source behaviour. It is also more likely 
to see a very fast flare in each source the longer they are moni- 
tored. More insight into the errors can be obtained when our new 
Dvar are compared to other studies (cf. Chapt. lTSI i. 

In Table|2]we show the median values of log(rb,var) and Dvar 
for each source class separately and for HPQs and LPQs com- 
bined together (FSRQ). The distributions are shown in Figs. |2] 
and[3]and we can see that the distributions of quasars and BLOs 
have considerable overlap. HPQs have a tail extending to higher 
Dvar and LPQs seem to be in between the HPQs and BLOs. Also, 
it is interesting to note that all the quasars are clearly Doppler- 
boosted, the smallest estimated Tb^var being 3.5 x 10" K for 
1928-1-738. We ran the Kruskal-Wallis analysis to examine the 
differences between source classes. (All Kruskal-Wallis analy- 
ses in this paper have been performed with the Unistat statistical 
package for Window^ (version 5.0).) The results confirm that 
all the source classes differ from the other classes significantly 
with a 95% confidence limit. 

3.2. Comparison witii previous analyses 

Our sample has 71 sources in common with the sample of LV99. 
We have re-calculated the Doppler boosting factors of LV99, us- 
ing the current cosmological model. Figure |4] shows the corre- 
lation between the Doppler boosting factors of LV99 and the 
new values. The results are very similar and the values correlate 
with a coefficient r - 0.77 (p=0.0000). Kruskal-Wallis analysis 
also shows that the values come from the same population. This 
confirms that the method is reliable because the results have not 
changed even though we have now ten more years of data. The 
differences are mainly due to poor fits in LV99 which are due 
to poor sampling or low flux density (large scatter) in the data. 
The scatter between the old and the new values is consistent with 
the error analysis in Sect. 13. II In some cases (e.g. 0430-t-052 and 
1 156+295) the source has clearly changed its behaviour and ex- 
hibits a much faster flare in our new dataset. The new estimates 
which are calculated using almost 30 years of data should there- 
fore be more representative of the source behaviour. 

We also compared our Dvar values with the Jorsta d et al.l 
( 120051) values for 15 AGN, obtained at 43 GHz. Figure |5] shows 
the correlation between the two values, and a Spearman rank cor- 
relation giv es a coefficient r=0.5 6 (p=0.0123). We can see that 
the values in lJorstad et al.l (l2005h are in general somewhat higher 
than ours. This can be due to their higher observing frequency 
(cf. Discussion). Also, their analysis method gives only an up- 
per limit to some sources. The distribution of source classes is 
similar to ours, with quasars having the highest Doppler boost- 
ing factors, GALs the lowest and BLOs being in between them. 



^ http://www.unistat.com/ 
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Fig. 2. Distribution of Tb.var of the fastest flare in each source 
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Fig. 3. Distribution of Dvar of the fastest flare in each source 



3.3. Core dominance 



Therefore we conclude that the Doppler boosting factors of these 
two analyses correspond well with each other. 

iHoman et alj ([2006) argued that during the most active state 
T'b.int should be closer to 2 x 10^'K and therefore the Dvar of the 
LV99 are overestimated. This would make our Dvar values even 
smaller, and the correspondence to Jorstad et al. ( 2005 ) would be 
worse. Higher rb,int would also increase our Lorentz factors for 
the fastest sources, and as is shown later (cf. Chapters]?] and |5]), 
our new values are akeady twice as high as in LV99 and in some 
sources even extremely high. We also note that a different value 
for Tb^int does not change the distributions themselves, only the 
numerical values. 

We also compared our Doppler boosting factors to a re - 
cent study at a lower frequency of 5 GHz (.Britzen et al.ll2007]) . 
They calculated the IC Doppler boosting factors by using VLBI 
data from the Caltech-Jodrell Bank Flat-Spectrum source sample 
jTaylor et al. 1996) and non-simultaneous ROSAT X-ray data. 
The Spearman rank correlation between our 24 common sources 
(excluding one outlier, 08364-710, with Die - 88) is still quite 
good (r=0.63, p=0.0004), and the slope of the linear fit between 
Die and Dvar is almost exactly one. We believe that most of the 
scatter is due to the errors in the IC Doppler boosting factors, 
since LV99 showed that for several reasons these are likely to be 
much less accurate than the variability Doppler boosting factors. 



Standard beaming models expect that more core-dominated ob- 
jects should be more beamed and thus have higher Doppler 




Fig. 4. Correlation between the Dvar from LV99 and our new Dvar 
values 
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Fig. 5. Correlation between the Dvar and (5var from iJorstad et alj 
(I2005h . 



boosting factors. We studied this by calculat ing the core- 
domin ance parameter R from VLBA data of iKovalev et alJ 
at 15 GHz. The core-dominance is calculated by relat- 
ing the flux density of the core ^ core to the total single-dis h flux 
density observed at 15 GHz 5 tot- In iKovalev et alj 32005b these 
are given for 250 sources observed with the VLBA at 15 GHz 
at different epochs. Their sample includes 80 sources for which 
we have determined Dvar- We calculated the median R for each 
source from the separate epochs. Figure |6] shows the correla- 
tion between logR and logDvai, excluding an outlier source 
0923+392 with a very small core dominance (log/? = -1.76). 
Spearman rank correlation between the parameters is r = 0.37 
(p = 0.0004). When the outlier source is included, the cor- 
relation is still significant with a coefficient r = 0.39 (p - 
0.0002). This shows that there is indeed some indication that 
so urces which ar e more core-dominated are also more boosted. 
In IKovalev et al.l (j2005) the core-dominance is defined to be the 
relation of S core to the total VLBA flux density S vlba- We tested 
the correlation using also this parameter but the results did not 
change because 5 tot and 5 vlba are so similar Usin g the core- 
dominance defined with S vlba, IKovalev et al.l (l2005i) show that 
quasars and BLOs are significantly different from GALs with 
lower core-dominance. 

Similar calculations were made in iBritzen et al.l ( |2007b for 
their sample. They used the core flux density and total single- 
dish flux density at 5 GHz to calculate the core-dominance pa- 
rameter They found no significant correlation between their 
core-dominance parameter and the IC Doppler boosting factor 
We have only 25 sources in common with their sample, and 
when we compared our Dvar with their core-dominance parame- 
ter, we found no correlation. 

4. The Lorentz factors and viewing angles 

We were able to calculate Lorentz factors Fvar and viewing an- 
gles ^var for 67 sources. The median values of different source 
classes are shown in Table [3] These are affected by two outliers 
with exceptionally large Fvar- The source 0923+392 has Fvar - 
216.1 and 1730-130 has Fvar = 64.6. Both of these show high 
superluminal motion of ySapp > 35c, which increases the Lorentz 
factors. The distributions of the source classes are shown in Figs. 
|7] and [8] In Fyar the distributions of quasars and BLOs over- 
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Fig. 6. Correlation between log(K) using data from lKovalev et all 
(l2005h and log(Z)var), excluding the outlier source 0923+392 
wifli \og(R) = -1.76 and log(Dvar) = 0.63. 

Table 3. Median values of Fvar and 0^..^. 



Type 


N 


r 

A var 




HPQ 


26 


17.41 


3.28 


LPQ 


23 


13.96 


3.90 




21° 


12.65 


3.96 


FSRQ 


49 


16.24 


3.37 


BLO 


13 


10.29 


5.24 


GAL 


5 


1.82 


15.52 


ALL 


67 


13.96 


3.81 



" = excluding outliers 0923+392 and 1730-130 



lap, but BLOs and GALs have slower jet speeds than quasars. 
Kruskal-Wallis analysis shows that when Fvar are studied with- 
out the outlier sources, HPQs differ from other classes with 
higher Lorentz factors, and GALs differ with smaller Lorentz 
factors. BLOs and LPQs come from the same population with 
a 13% confidence. There is one BLO (1823+568) for which 
Fvar = 37.8. This source has a lso been classified as a HPQ 
(e.g. lVeron-Cettv & Ver"onll2006l) and therefore we ran the KW- 
analysis again by moving this source into the HPQ class. In this 
case also BLOs and LPQs differ significantly from each other 
with a 96% confidence. We must note that in our samples of 13 
to 26 objects the significance of differences can depend on the 
classification of a single extreme source. However, a clear re- 
sult is that the BLOs and the quasars differ from each other with 
BLOs having slower jets (FSRQs and BLOs differ significantly 
with a 99% confidence if 1823+568 is classified as a BLO and 
with a 99.9% confidence if 1823+568 is classified as a HPQ). 
This result is in accordance with several earlier but more indirect 
estimates of jet s peeds. Similar resul ts have also been obtained 
with simulations dHughes et al.ll2002l) . 

When 0var is studied, the distributions overlap even more 
and KW-analysis shows that GALs and BLOs differ from other 
source classes with a 95% confidence. Also, if the differences be- 
tween FSRQs and BLOs are studied, they differ from each other 
significantly with a 99% confidence with BLOs having larger 
viewing angles. 
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Fig. 7. Distribution of Fvai of the fastest flare in each source (ex- 
cluding the outlier LPQs 0923+392 and 1730-130) 



We compared our new Lorentz factors with the ones from 
LV99. There are 38 sources in common in our samples, and in 
Fig. |9] we can see that the new Fyar are about twice as large as 
in LV99. This is mainly because also the apparent speeds used 
in our analysis are about twice as large as in LV99, where they 
were collected from the literature. We should now have better 
and more homogeneous estimates of ySapp which should make 
our new estimates more accurate. The correlation between our 
Fvar and those from LV99 is only r = 0.37 (p - 0.0113) when 
the oudier (0923+392) is not included (the other outlier 1730- 
130 is not included in the sample of LV99 ). When we c o mpare 
our Fvar with Lorentz factors estimated in iJorstad et al.l (l2005h 
we find a good correlation of r = 0.50 (p = 0.0281) which gives 
credibility to our new estimates. Also the distribution is similar 
with quasars having the highest Lorentz factors, BLOs being in 
between and GALs having the smallest Lorentz fac tors. We have 
only 1 2 sources in common with the sample of iBritzen et al.l 
(l200 7l) and we find no correlation between our Lorentz factors 
even if we leave out their outlier source 0016+73 1 for which they 
determine a Lorentz factor of 860. The differences are probably 
due to their much lower values of ySapp and differences in the 
Doppler boosting factors. 

When studying the viewing angles, the difference between 
our new and LV99 values is not as large and the correlation is 
goo d (r = 0.51 p — 0.00 05). The correlation between our values 
and lJorstad et al.l (l2005i) values is also very good with r - 0.59 



Fig. 8. Distribution of ^var of the fastest flare in each source 



(p - 0.0105). Again ther e are no significant correlations be- 
tween Avar and values from lBritzen et al.l (|2007[) . 



5. Discussion 

Although the Doppler boosting factors have remained on the 
average almost identical (cf. Fig. |4| even though we now have 
ten more years of data, there is a factor of two difference in the 
Lorentz factors when comparing our new results to LV99. This is 
mainly due to twice as fast apparent speed s in our new analysis. 
In LV 99 the ySapp values are mainly from IVermeulen & CohenI 
(Il994h . in which the values were gathered from the literature. 
This means that the speeds were obtained in the 1980s and early 
1990s from a variety of observing programmes and frequen- 
cies, mainly at 5 GHz. We use 15 GHz MOJAVE data from 1994 
up to September 2008, giving a far more uniform dataset. We 
also note that higher frequency VLBI observations often tend to 
give high er apparent speeds as can b e seen by comparing, e.g ., 
data from Vermeul en & Cohe n (1994'):'Kenerm ann et al](l2004l) : 
Jorstad et al. (20051) and^itzen et al. (2007), although the rea- 
son for this is not clear. The use of a single frequency therefore 
diminishes the internal scatter. 

Long-term observations, both TFD monitoring and VLBI, 
are essential in understanding source behaviour. We have al- 
ready shown this in our earl ier studies of long-term variability 
behavi our of these sources dHovatta et al.l 120071: iHovatta et al.l 
l2008allbl) . Some sources have changed their behaviour in our 
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Fig. 9. Correlation between Lorentz factors from LV99 and 
the new values (excluding the outlier source 0923+392 with a 
Lorentz factor of 216 in the new analysis). 



TFD observations during the past 10 years. Similar changes can 
be seen in VLBI sources over a lon g time, such as the pos ition 
swings seen, for example, in 3C 273 (ISavolainen et alj2006h and 
NRAO 150 dAgudo et al.ll2007h . We emphasise that our way to 
estimate the Lorentz factors and viewing angles is as good as can 
be, when we are characterising a complex, changing jet with just 
two parameters, both assumed to be constant. 

In addition to having twice as high apparent speeds in gen- 
eral, there are also some extreme sources showing very high ap- 
parent motion. While in LV99 the highest apparent speed was 
14.9, we now have two sources (0923+392 and 1730-130) with 
jSapp > 35c. As a consequence, the Lorentz factors of these 
sources are also extremely high, Fvar - 216 for 0923+392 and 
Tvar = 65 for 1730-130. We also note that if the yS^pp = 45.9 for 
the source 1510-089 from Jorstad et al. (2005) is accepted, our 
log(rb,var) - 14.4 would indicate a Fvar =71. Thus, the existence 
of a class of very fast jets should perhaps at least be considered 
as a possibility. 

However, at least the Fvar = 216 seems rather unlikely in 
view of our current knowledge of the jets in AGN. One alterna- 
tive explanation is that 0923+392 has a higher observed bright- 
ness temperature than what we have obtained (log(rb,var) = 12.6) 
from our monitoring. If we saw changes of about 1 Jy within a 
time period of a week, the brightness temperature would be of 
the order of lO'^/T, which would change the Lorentz factor to a 
more acceptable value of under 50. Our sampling, like in other 
monitoring programmes, is too sparse to detect such rapid flares 
reliably, and therefore we have initiated a denser monitoring 
schedule for the source. Another possibility is that the source has 
a complicated internal structure or geometry, such that the TFD 
variations an d the apparent speeds do not refer to t he same com- 
ponent (e.g. l Alberdi et alJll993l: iFev et alJll997t lAlberdi et alj 
[2000). However, independent of our Tb^var estimates, 0923+392 
must have a Lorentz factor of at least 43 as F > ySapp. 

In Fig. [To] we have plotted the observable quantities ySapp 
and log(rb,var)- We have included curves to mark areas of dif- 
ferent Fvar and Avar- The outlier sources are clearly visible in 
this plot but otherwise the sources are within rather well-defined 
limits. Almost all the sources have Fvar < 40 and 0var < 20°. 
The differences between the source classes are also seen in this 
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Fig. 10. Observable quantities ySapp and log(rb,var) together with 
intrinsic parameters Fvar and ^var- 



plot, with GALs having slow speeds and low brightne s s tem- 
peratures. Using Monte-Carlo simulations, ICohen et al.l ( l2007h 
find an upper limit for the Lorentz factor to be F « 32. This 
agrees quite well with our results, although in our sample we 
have five sources with Fvar between 30 and 50 and two quasars 

with Fvar > 50. 

A common assumption is that sources are viewed close to 
the critical angle 9c = 1/F. Figure [TT] shows the Fsin0 distribu- 
tion for our sources, which indeed does peak around 1 . However, 
a number of sources have F sin 6 significantly larger or smaller, 
so that the assumption of F sin = 1 will in many cases lead 
to false conc l usions both for individual sources and for samples. 
ICohen et alj (l2007h have calculated a Monte Carlo simulation 
for a flux de nsity limited survey (similar to ours), (Fig. Ic in 
ICohen et all (2007), in which the sources with F = 15 are plot- 
ted), finding a distribution rather similar to our empi rical one. 

The simplest unific ation scheme for AGN (e.g. B arthelll989l : 
lUrrv & Padovanill995h predicts that in all radio quasars (includ- 
ing ordinary quasars and FSRQs) the F distributions should be 
similar. Ordinary quasars, on the other hand, should have larger 
viewing angles than blazars (FSRQs and BLOs). BLOs could 
also have a different F distribution because their parent popula- 
tion is different than in quasars. In Fig. [12] we show the Lorentz 
factors and viewing angles in a polar plot. It is easy to see that al- 
most all the sources are seen at a small viewing angle. What must 
be kept in mind is that our sample essentially includes the ~ 100 
brightest northern compact radio sources, so this is as expected 
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Fig. 11. Distribution of Fsinfl, excluding the source 0923+392 

with Tvar sin 0var =10. 



for Doppler boosting dominated sources. Although we have di- 
vided our quasars into HPQs and LPQs, they all are FSRQs (or 
blazars), and we do not see ordinary quasars with 0var > 20°. 
Similarly, all our BLOs are radio se lected BLOs (RBLs) which 
should have small viewing angles (l Barthel|[T98 9). In addition, 
the five GALs in our sample are the most compact and variable 
sources of their type. Figure [13] also shows how the sources at 
different redshifts are selected due to their Doppler boosting so 
that small Doppler boosting factors are not seen at high redshifts. 
These caveats must be kept in mind when comparing our differ- 
ent classes of sources with each other. 

The d istribution of Fyar agrees well with unification sc hemes 
for AGN. iPadovani & Uml j 19921) and lUrrvet all (|1991|) calcu- 
lated beamed luminosity functions for FSRQs and BLOs, and 
compared them with observed luminosity functions, which were 
assumed to have the same shape as the intrinsic luminosity func- 
tion of radio-loud quasars. The validity of this assumption was 
recently confirmed by using the maximum l ikelihood method to 
calculate the intrinsi c lu minosity functions dLiu & Zhang||2007l; 
ICara & Listed |2008|) . In |Padovani &_0n3 il99^ the FSRQs 
were best described with a distribution of 5 < F < 40 with an 
average of ~ 1 1 (using a cosmology with Hq = 50 km s"' Mpc"' 
and qo = 0). In our sample all the FSRQs have Fvai > 5 and three 



have Fvar > 40. The median Fyar is som ewhat bigg e r than the 
average from the luminosity functions. In Urrv et all (|1991|) the 
BLOs were best described with a distribution of 5 < F < 35 with 
an average of ~ 7. Again we have only three BLOs with Fyai < 5 
and one BLO(1823-h568, also classified as HPQ) with Fyar > 35. 
Our median is also only slightly bigger. In addition they made 
another fit for the BLOs with a distribution of 2 < F < 20 which 
also fit the data well. In this case only one BLO (1807+698 with 
Fyar = 1 .0) would be below the lower limit and there would be 
two BLOs with Fyar > 20. Their calculations also set a value 
for the critical viewing angle in which all the FSRQs and BLOs 
should be seen. For FSRQs the angle is 14°, and all our quasars 
are within this limit. For BLOs the critical angle is either 11° 
or 19° depending on the model. In both cases only one BLO 
(1807+698 with 6lyai = 57°, which is also the nearest BLO in our 
sample) would have a viewing angle outside the limit. 

As expected with sources with small viewing angles, distri- 
butions of LPQs and HPQs overlap and one should not iden- 
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Fig. 12. Polar plot of Lorentz factors against the viewing an- 
gles, excluding the outlier sources 0923+392 and 1730-130. 
The lower panel shows only the 0° to 20° portion of the view- 
ing angles, additionally excluding the galaxies 0007+106 and 
0316+413, and the BLO 1807+698. 



tify LPQs as ordinary quasars. In Fig. [14] we show the maxi- 
mum optical polarisation against the viewing angle for all the 
sources for which we could find polarisation information in the 
literature (f max and its reference for all the sources are shown 
in Table [TJ. Most of the sources seem to be within an enve- 
lope with Pmax decreasing as 0yai- increases. The situation may 
be similar as in Se yfert galaxies with an obscuring torus (e.g. 
ISchmittetalJl2001h . When the viewing angle is small, we can 
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Fig. 14. Maximum optical polarisation against the viewing an- 
gle. 



see deeper into the jet and are more likely to see the source 
as more optically polarised. Also, it has been shown that both 
the optical and the radio polarisation originate from transverse 
shocks located very close to the base of the jet (iLister & SmithI 
12000 ). Therefore the range of polarisation variability is much 
higher for sources with a viewing angle of a few degrees than 
for sources with 0var = 10° - 20°. In Fig. [14] it is also pos- 
sible to see that the definition P^ax = 3% for an object to be 
classified as highly polarised is somewhat artificial. We believe 
that with a larger number of polarisation observations, a con- 
siderable fraction of the LPQ population could, at times, show 
polarisation exceeding 3%, as expected from blazar-type AGN. 
Since there usually is only a small number of polarisation obser- 
vations, the likelihood of observing high polarisation increases 
wi th decreasing viewing angle. However, as has been pointed out 
bv lLister & SmithI (120001) . there is a physical difference between 
at least some HPQs and LPQs, especially in the magnetic field 
structure of their inner jets. 

Figure [14] hints that the transition from blazar-type AGN to 
ordinary quasars may occur around a viewing angle of 15° to 
20°, presumably corresponding to the half opening angle of an 



AGN obscuring torus, and in agreement with, e.g., estimates 
from luminosity functions and source counts. As our sample 
does not contain any ordinary quasars, no definite conclusions 
can be drawn. 



6. Conclusions 

We have decomposed flux curves of 87 sources into exponential 
flares and using the fits calculated the variability brightness tem- 
perature and Doppler boosting factor for each source. In addi- 
tion we used new MOJAVE observations of apparent jet speeds 
to calculate the variability Lorentz factors and viewing angles. 
We have compared our results with LV99, in which the param- 
eters were determined in a similar way. In our new analyses we 
have used almost 15 years more data and the estimates should 
be more accurate. Our main conclusions can be summarised as 
follows. 

1 . The variability Doppler boosting factors have remained on 
the average almost identical compared to LV99. All the 
quasars are Doppler-boosted and they are in general more 
boosted than the BLOs or GALs. 

2. The Lorentz factors of our new analyses are about twice as 
large as in LV99. The difference can be explained with twice 
as large apparent VLBI speeds in our new analyses. The 
BLOs have slower jets compared to quasars (Fblo = 10.3, 
Tfsrq = 16.2 and Fgal = 1-8). 

3. A few sources, 0923-1-392 and 1730-130, have extreme 
Lorentz factors (216 and 65, respectively). Either these 
Lorentz factors are real, or the sources exhibit so rapid flares 
that the fast variations have remained undetected in the mon- 
itoring programmes. We are studying the second possibility 
with dense observations of 0923+392. A third possibility is 
that the sources have a structure too complex for our method. 

4. Almost all the sources in our sample are seen in a small view- 
ing angle of less than 20 degrees (0blo = 5.2°, 0fsrq = 3.4° 
and 0GAL = 15.5°). 

5. The viewing angle distribution peaks around Fsinfl = 1, 
with a distribution similar to that found in simulations. 

6. Our results generally follow the predictions of basic unifi- 
cation models for AGN. Based on our results, we cannot 
separate HPQs and LPQs from each other, and therefore it 
is well-grounded to treat them as a single group of FSRQs 
when the jet parameters are considered. 
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Table 1. continued. 



B1950-name 


Other name 


Type 


V 


z 


Quality 


logCn) 




j8a 


Evai 




logR 




ref. 


1823+568 


4C 56.27 


BLO 


22 


0.663 


A 


13.11 


6.4 


20.950 


37.8 


5.0 


-0.154 


29.8 


5 


1828+487 


3C 380 


LPQ 


37 


0.692 


A 


12.97 


5.7 


13.66 


19.3 


7.1 


-0.391 


0.4 


2 


1928+738 


4C 73.18 


LPQ 


22 


0.302 


A 


11.55 


1.9 


8.48 


19.9 


12.8 


-0.402 


0.8 


2 


1954+513 




LPQ 


22 


1.223 


A 


13.30 


7.4 








-0.181 


1.4 


11 


2005+403 




QSO 


37 


1.736 


A 


14.07 


13.3 


19.48 


21.0 


4.0 


-0.257 






2007+776 


S5 2007+77 


BLO 


22 


0.342 


A 


13.39 


7.9 








-0.304 


15.1 


11 


2021+614 


OW637 


GAL 


22 


0.227 


A 


11.65 


2.1 


0.41 


1.3 


13.3 




0.3 


3 


2022+171 




LPQ 


22 


1.050 


A 


12.99 


5.8 














2121+053 




HPQ 


22 


1.941 


A 


14.25 


15.3 


13.05 


13.2 


3.7 


-0.061 


10.7 


7 


2134+004 


OX 057 


LPQ 


22 


1.932 


A 


14.32 


16.1 


5.62 


9.0 


2.2 


-0.405 


2.8 


12 


2136+141 




LPQ 


22 


2.427 


A 


13.46 


8.3 


5.08 


5.8 


6.2 


-0.175 


1.6 


11 


2145+067 




LPQ 


37 


0.990 


G 


14.28 


15.6 


2.49 


8.0 


1.1 


-0.171 


0.6 


2 


2200+420 


BLLAC 


BLO 


37 


0.069 


E 


13.28 


7.3 


10.700 


11.6 


7.3 


-0.264 


23.0 


1 


2201+315 


4C 31.63 


LPQ 


22 


0.295 


G 


13.17 


6.7 


7.88 


8.1 


8.5 


-0.096 


0.2 


2 


2223-052 


3C446 


HPQ 


22 


1.404 


G 


14.31 


16.0 


16.470 


16.5 


3.6 


-0.203 


17.3 


10 


2227-088 




HPQ 


22 


1.562 


A 


14.31 


15.9 


4.95 


8.8 


2.0 


-0.137 


6.9 


12 


2230+114 


CTA 102 


HPQ 


37 


1.037 


E 


14.28 


15.6 


15.51 


15.5 


3.7 


-0.269 


10.9 


7 


2234+282 




HPQ 


22 


0.795 


A 


13.03 


6.0 








-0.464 


4.4 


11 


2251+158 


3C 454.3 


HPQ 


37 


0.859 


E 


15.26 


33.2 


14.86 


19.9 


1.3 


-0.564 


16.0 


1 


2254+074 


PKS 2254+074 


BLO 


22 


0.190 


A 


11.91 


2.5 










21.0 


1 



Refere nces: (1)1 Angel & Stockmanl(ll980l): qillmpey & TapiaHl990h: G i llniDev et al.Hl991h; (4) llmpev et"aLH2000h; ( 5) lJorstad et al](l2007D: 
(e'jiLister & Smith (2000); f7) Moore & Stockmanl i ll984l) : (8) lPunsH(ll996l) : ('9) ISluse et al.l(l2005h : ri0) lsmith et alJ jl985h : (ll) IStickel & KuehJ 
jl994l) : (12) .Willsetal., H992i) . 
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